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We use 2 H NMR stimulated-echo spectroscopy to measure two-time correlation functions charac- 
terizing the polymer segmental motion in polymer electrolytes PPO-LICIO4 near the glass transition 
temperature T g . To investigate effects of the salt on the polymer dynamics, we compare results for 
different ether oxygen to lithium ratios, namely, 6:1, 15:1, 30:1 and oo. For all compositions, we 
find nonexponential correlation functions, which can be described by a Kohlrausch function. The 
mean correlation times show quantitatively that an increase of the salt concentration results in a 
strong slowing down of the segmental motion. Consistently, for the high 6:1 salt concentration, a 
high apparent activation energy _E a = 4.1 eV characterizes the temperature dependence of the mean 
correlation times at T g < T < 1.1 T g , while smaller values E a ~ 2.5 eV are observed for moderate 
salt contents. The correlation functions are most nonexponential for 15:1 PPO-LiClCU, whereas 
the stretching is reduced for higher and lower salt concentrations. This finding implies that the 
local environments of the polymer segments are most diverse for intermediate salt contents, and, 
hence, the spatial distribution of the salt is most heterogeneous. To study the mechanisms of the 
segmental reorientation, we exploit that the angular resolution of 2 H NMR stimulated-echo exper- 
iments depends on the length of the evolution time t p . A similar dependence of the correlation 
functions on the value of t p in the presence and in the absence of ions indicates that addition of salt 
hardly affects the reorientational mechanism. For all compositions, mean jump angles of about 15° 
characterize the segmental reorientation. In addition, comparison of results from 2 H and 7 Li NMR 
stimulated-echo experiments suggests a coupling of ion and polymer dynamics in 15:1 PPO-LiClCU. 



I. INTRODUCTION 

The glass transition of polymer melts is of tremendous 
interest from the viewpoints of fundamental and applied 
science. This importance triggered a huge amount of 
studies, which established that the glass transition phe- 
nomenon of polymers is related to the a processi 2 ^ An 
improvement of our understanding of the glass transition 
thus requires detailed characterization of the microscopic 
dynamics underlying the a process. In particular, it is 
important to determine both the time scale and the ge- 
ometry of the molecular dynamics. 

Various methods were used to measure the correlation 
time of the a process, r Q . For polymers, a striking feature 
of the a process is a non-Arrhenius temperature depen- 
dence. In most cases, the Vogel-Fulcher-Tammann equa- 
tion, r a (T) = t exp [B/(T — To)), enables a satisfactory 
interpolation of the data, suggesting a divergence of the 
correlation time at a temperature To not too far below 
the glass transition temperature T g , i.e., To«T g — 50 K. 3 
Hence, high apparent activation energies E a character- 
ize the temperature dependence of the a process close 
to T y . xi Another intriguing property of the a process is a 
nonexponential loss of correlation, which is often well de- 
scribed by the Kohlrausch- Williams- Watts (KWW) func- 
tion, exp[— (t/r)P]. By contrast, few experimental tech- 
niques yield detailed insights into the geometry of the 
molecular motion associated with the a process. 

2D NMR provides straightforward access to both the 
time scale and the geometry of molecular reorientation 
in polymeric and molecular glass-forming liquids. 3 ' 4 ' 5,6 '- 7 
2D NMR experiments can be performed in the frequency 



domain 8 ! 9 ! 10 ! 11 and in the time domaim 12 i 13 i 14 i 15 i 16 i 17 i 18 
The latter approach is also called stimulated-echo (SE) 
spectroscopy. Concerning the geometry of polymer dy- 
namics, 2D NMR studies ruled out the existence of well- 
defined jump angles 7 for the rotational jumps associ- 
ated with the a processi 3 ^LSi24i2i 1 ii 8 . Instead, results 
from 2D and 3D NMR showed that the reorientation is 
comprised of rotational jumps about angles 7 = 10 — 40° 
between which the polymer segments show isotropic ro- 
tational diffusion, 11 ! 18 ! 19 ! 20 ! 21 

We carry out 2 H NMR SE experiments to investigate 
the segmental motion of poly(propylene oxide) (PPO) 
containing various amounts of the salt LiC104. Thus, we 
focus on polymer electrolytes, which are promising candi- 
dates for applications in energy technologies. 22 Since the 
use of polymer electrolytes is still limited by the achiev- 
able electric conductivities, there is considerable interest 
to speed up the ionic diffusion. It has been established 
that segmental motion of the host polymer facilitates the 
ionic migration^ 2 ^* 2 * yet the understanding of this cou- 
pling is incomplete. In dielectric spectroscopy, the ionic 
diffusion coefficient and the relaxation rate of polymer 
segments in contact with salt were found to show a linear 
relationship,- 25 suggesting that segmental motion triggers 
the elementary steps of the charge transport. 

Polymer electrolytes PPO-LiC104 have attracted much 
attention because of high ionic conductivities together 
with their amorphous nature^ When studied as a 
function of the salt concentration, the electric conduc- 
tivity Ode at room temperature shows a broad max- 
imum at ether oxygen to lithium ratios, 0:Li, rang- 
ing from about 10:1 to 30:lm^2& For such intermedi- 
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ate salt concentrations, differential scanning calorimetry 
showed the existence of two glass transition steps 
implying a liquid-liquid phase separation into salt-rich 
and salt-depleted regions. An inhomogeneous structure 
was also assumed to explain results from dielectric and 
photon correlation spectroscopies for intermediate salt 
concentrations, 2 ^' 3 - which give evidence for the pres- 
ence of two processes related to fast and slow poly- 
mer dynamics. By contrast, the existence of salt-rich 
and salt-depleted microphases was challenged in neutron 
scattering work on mixtures PPO-LiClC^iSi*^ 2 " Finally, 
it was proposed that a structure model of 16:1 PPO- 
LiClC>4 resulting from reverse Monte-Carlo simulations 
of diffraction data enables reconciliation of the apparent 
discrepancies^ It does feature salt-rich and salt-depleted 
regions, but the structural heterogeneities occur on a 
small length scale of 1 nm. 

In NMR work on polymer electrolytes, the ion dynam- 
ics was studied in some detailj 3 ^ 3 ^ 3 ^^ 3 ^ 3 ^* 3 ^ 4 ^ e.g., ionic 
self diffusion coefficients were determined using field- 
gradient techniques while investigations of the 
polymer dynamics were restricted to analysis of the spin- 
lattice relaxation. 38 ) 40 ! 43 Only recently, modern NMR 
methods were employed to study the segmental motion 
in organic-inorganic hybrid electrolytes^ We use vari- 
ous 2 H and 7 Li NMR techniques for a detailed charac- 
terization of the polymer and the ion dynamics in poly- 
mer electrolytes PPO-LiC104. Our goals are twofold. 
From the viewpoint of fundamental science, we ascertain 
how the presence of salt affects the polymer dynamics. 
In this way, we tackle the question how structural het- 
erogeneities introduced by the addition of salt affect the 
glass transition of polymers. From the viewpoint of ap- 
plied science, we investigate how the polymer assists the 
ionic diffusion so as to give support to the development 
of new materials. For example, we compare 2 H and 7 Li 
NMR data to study the coupling of ion and polymer dy- 
namics. 

Very recently, we presented results from 2 H and 7 Li 
NMR line-shape analysis for polymer electrolytes PPO- 
LiC10 4 i£ The 2 H NMR spectra indicate that addition of 
salt slows down the segmental motion. Moreover, the line 
shape shows that broad distributions of correlation times 
G(log r) govern the polymer dynamics. These dynamical 
heterogeneities are particularly prominent for intermedi- 
ate salt concentrations. While such line-shape analysis 
provides valuable insights on a qualitative level, the pos- 
sibilities for a quantification are limited. 

Here, we exploit that 2 H NMR SE spectroscopy en- 
ables measurement of two-time correlation functions that 
provide quantitative information about polymer segmen- 
tal motion on a microscopic levels In this way, we con- 
tinue to investigate to which extent presence of salt af- 
fects the a process of polymer melts. While previous 
2 H NMR SE studies focussed on molecules featuring a 
single deuteron species, two deuteron species, which ex- 
hibit diverse spin-lattice relaxation (SLR) and spin-spin 
relaxtion (SSR), contribute to the correlation functions 



in our case. Therefore, it is necessary to ascertain the in- 
formation content of this method in such situation. We 
will demonstrate that, in the presence of two deuteron 
species, it is necessary to use appropriate experimental 
parameters, when intending to correct the measured data 
for relaxation effects, which is highly desirable in order 
to extend the time window of the experiment. 

II. THEORY 

A. Basics of 2 H and 7 Li NMR 

Solid-state 2 H NMR probes the quadrupolar preces- 
sion frequency luq, which is determined by the orienta- 
tion of the quadrupolar coupling tensor with respect to 
the external static magnetic flux density Bo- This ten- 
sor describes the interaction of the nuclear quadrupole 
moment with the electric field gradient at the nuclear 
site. The monomeric unit of the studied PPO, [CD2- 
CD(CD3)-0]„, features three deuterons in the backbone 
and three deuterons in the methyl group, which will be 
denoted as B deuterons and M deuterons, respectively. 
For both deuteron species, the coupling tensor is approx- 
imately axially symmetric^ Then, u>q is given by 3 * 

u Q (9) = ±^(3cos 2 9 y -l) (1) 

Here, y = B,M refers to the different deuteron species. 
For the B deuterons, 9b specifies the angle between the 
axis of the C-D bond and Bo and the anisotropy parame- 
ter amounts to <5s«27rxl20kHz.— For the M deuterons, 
the coupling tensor is preaveraged due to fast threefold 
jumps of the methyl group at the studied temperatures. 
Consequently, 9m is the angle between the threefold axis 
and Bo and the anisotropy parameter is reduced by a 
factor of three, 8m — <5b/3i£ The two signs in Eq. (JJJ 
correspond to the two allowed transitions between the 
three Zeeman levels of the 1—1 nucleus. We see that, for 
both deuteron species, the quadrupolar precession fre- 
quency is intimately linked to the orientation of well de- 
fined structural units and, hence, segmental motion ren- 
ders ojq time dependent, which is the basis of 2 H NMR 
SE spectroscopy^ 

In solid-state 7 Li NMR (7 = 3/2), the quadrupolar in- 
teraction affects the frequencies of the satellite transitions 
±3/2 <-» ±1/2, while that of the central transition 1/2 <-> 
— 1/2 is unchanged. For polymer electrolytes, diverse 
lithium ionic environments lead to a variety of shapes of 
the quadrupolar coupling tensor so that 7 Li NMR spectra 
exhibit an unstructured satellite component Analy- 
sis of 7 Li NMR SE experiments will be straightforward if 
the time dependence wq(<) is solely due to lithium ionic 
motion. This condition is approximately met in appli- 
cations on solid-state electrolytes, i.e., when the lithium 
ionic diffusion occurs in a static matrix^&ii For poly- 
mer electrolytes, the lithium environments change due 
to ion and polymer dynamics, occurring on comparable 
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time scales. Consequently, both these dynamic processes 
can render the electric field gradient at the nuclear site 
and, hence, ojq time dependent^ Thus, when 7 Li NMR 
is used to investigate charge transport in polymer elec- 
trolytes, it is a priori not clear to which extent the lithium 
ionic diffusion and the rearrangement of the neighboring 
polymer chains are probed. 



B. 2 H NMR stimulated-echo spectroscopy 

2 H NMR SE spectroscopy proved a powerful tool to 
study molecular dynamics with correlation times in the 
range of 10 fis < r < 1 3,6,12,13,14,15,16,17,49,50,51,52,53 In this 

experiment, the frequency u)q is probed twice during two 
short evolution times t p -c t, which embed a long mixing 
time t m ps r. In detail, appropriate three-pulse sequences 
are used to generate a stimulated echoi 3 i 4 i 5 i 6 i 7 Evaluating 
the height of this echo for various mixing times t m and 
constant evolution time t p , it is possible to measure two- 
time correlation functions. Neglecting relaxation effects, 
one obtains depending on the pulse lengths and phases 
in the three-pulse sequenc o 3 ' 4 ' 5 ' 6 ' 7 

E^ s {t m ;t p ) = (sin[w Q (0)tp] • sm[u Q (t m )t p ] ), (2) 



p] ■ cos[u>Q(t m )t p ] ). (3) 

Here, the brackets (. . .} denote the ensemble average, in- 
cluding the powder average for our isotropic distribution 
of molecular orientations^ We see that the frequencies 
u>q(0) and u)Q(t m ) are correlated via sine or cosine func- 
tions. In general, molecular reorientation during the mix- 
ing time leads to ojq(0) ^ujQ(t m ) and, hence, to decays 
of E% (t m ) and E% (t m ). Normalizing these correlation 
functions, we obtain (x = SS,CC) 



Ei {tm j tp] 



E% {tra'i tp) 

EWvtp) 



(4) 



to which we will refer as sin-sin and cos-cos correlation 
functions in the following. 

The dependence of F 2 (t m ;t p ) on the evolution time 
provides access to the jump angles of molecular rotational 
jumps i 15 i 16 ' 17 i 18 ' 51 i 52 ' 53 This approach exploits that the 
angular resolution of the experiment depends on the 
length of the evolution time t p . When the quadrupolar 
coupling tensor is axially symmetric, the sin-sin correla- 
tion function for t p — > yields the rotational correlation 
function of the second Legendrc polynomial, F 2 (t m ): 

E 2 {tm , tp ->-0) cx tl(uj Q (0)ujQ(t m ) ) (X 

cx (P 2 [cos(0(O)]P 2 [cos(0(t m )]) oc F 2 {t m ). (5) 

On extension of the evolution time, smaller and smaller 
changes of the resonance frequency, |u>g(i m )— u>q(0)\, are 
sufficient to observe a complete decay of F 2 (t m ; t p ), i.e., 



isotropic redistribution of the molecular orientations is 
no longer required. In the limit t p — >oo, 



F 2 x (t 



d) = (S[w Q (t m )-u Q (0))) = F 2 J (t m ) (6) 



and, hence, the jump correlation function F 2 (t m ) is mea- 
sured, which characterizes the elementary steps of the 
reorientation process. The mean jump angle 7 of the 
molecular rotational jumps can be determined by com- 
parison of the correlation times tc and rj characterizing 
the decays of F 2 (t m ) and F 2 (t m ), respectively^ 



t.j *J . 2 
— = o sm 7- 

T C I 



(7) 



Analysis of the complete dependence of the correlation 
time on the evolution time, T(t p ), provides further infor- 
mation about the distribution of jump angles. 1617,18 

Here, the existence of two deuteron species complicates 
the analysis of sin-sin and cos-cos correlation functions. 
Thus, prior to a presentation of experimental results, it is 
important to determine the interpretation of these corre- 
lation functions in such situation. In the first step, we will 
continue to neglect relaxation effects and study conse- 
quences of the different anisotropy parameters 5b = 3<5m- 
In the second step, we will investigate how the interpre- 
tation of the correlation functions is affected by the fact 
that the deuteron species exhibit different spin-lattice 
and spin-spin relaxation times, T\ and T 2 , respectively. 

In the first step, the non-normalized correlation func- 
tions can be written as the following sum of contributions 
from the B and the M deuterons 



E 2 {t m ,tp) S TUyB^ ,y^tm'i tp) 



'2, y 



(t m ;t p ) = [a x y (t p )-z y '(t p )]<f ) *(t m ;tp) + z*(t p ) (8) 



Here, we consider that e 2 y (t m ]t p ) decays from an initial 
correlation aj(i p ) to a residual correlation Zy(t p ) and, 
hence, the correlation functions (j)y{t m \t p ) are normal- 
ized according to <j>y{0) = 1 and </>^(oo) = 0. The molar 
ratios of the deuteron species, m y = 1/2, are proportional 
to their contributions to the equilibrium magnetization. 
The initial state values of the correlation functions are 
determined by (£ = sin, cos) 



\, y (0;t p ) = al(t p ) = (e[u Q tp]) % 



(9) 



where (. . .) y denotes the powder average for the respec- 
tive species. The final state values depend on the ge- 
ometry of the reorientation processi 4 ^^ For an isotropic 
motion, the frequencies wq(0) and ujQ.(t m ) are uncorre- 
cted in the limit t m 3>T, leading to 



4, y {oo;t p ) = z x y (t p ) = {£[uj Q t p })l 



(10) 



In Fig. ^ we show theoretical values of the initial and 
residual correlations a^(t p ) and Zy(t p ), respectively. For 
the calculation, we used typical anisotropy parameters 
S B = 2tt x 120 kHz and 5 M = 2tt x 40 kHz. 45 Due to the 
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FIG. 1: Initial correlations a,y{t p ) and residual correlations 
Zy(t p ) for an isotropic motion, see Eq. ®. The values were 
calculated for anisotropy parameters 8b — 2n x 120 kHz and 
<5 M = 27rx4Q kHz. 

periodicity of the sine and cosine functions, all the quan- 
tities show an oscillatory behavior with a frequency pro- 
portional to the anisotropy parameter. Thus, as a conse- 
quence of Sbj^Sm, the ratios a x B /a x M are functions of the 
evolution time, indicating that the relative contributions 
of the two deuteron species to the correlation functions 
depend on the value of t p . 

SLR and SSR limit the lengths of the interpulse delays 
in SE experiments. To consider relaxation effects in the 
second step, we modify Eq. JHJ according to 

^2^m', tp) =2_j m y e 2,y(tm', tp)Rl y y(t m )R2 y y(2t p ) (11) 
V 

While Ri, y (t m ) describes the decrease of the echo height 
due to SLR during the mixing time, R2, y (2t p ) quantifies 
the loss of signal due to SSR during the evolution times, 
i.e., Ri iy (0) — 1 and Ri iV (oo) = (i = 1,2). In analogy 
with Eq. QJ, we define T% (t m ;t p ) = £2 {t m ]tp)/£i (0;i p ). 
The symbols £f and T% are used to stress that the cor- 
relation functions are affected by relaxation effects. As 
a consequence of fast methyl group rotation at the stud- 
ied temperatures, the SLR times of the M deuterons are 
shorter than that of the B deuterons, Tim ^ T\,b*^ 
Furthermore, the deuteron species exhibit different SSR 
and, hence, the ratio of the respective contributions to 
the correlation functions increasingly deviates from the 
theoretical value a x B (t p ) / a x M (t p ), when the evolution time 
is increased, see Appendix. 

The time window of the experiment is limited in par- 
ticular by the fast SLR of the M deuterons at the studied 



temperatures, T^m^ 10 ms. Thus, correction for SLR is 
desirable so as to extend the time window. In 2 H NMR, 
such correction is not possible for the sin-sin correlation 
functions, since the decay of Alignment order, which is 
present during the mixing time of this experiment, is not 
accessible in an independent measurement so that the 
damping due to SLR is unknown. By contrast, for the 
cos-cos correlation functions, the decay of Zeeman order, 
which exists during the mixing time of this measurement, 
can be determined in a regular SLR experiment. In pre- 
vious 2 H NMR SE studies, it thus proved useful to focus 
on the cos-cos correlation function and to divide by 
the SLR function Ri(t m ) to remove SLR effects. 

Following an analogous approach, we analyze the cor- 
relation functions 

£ cc (t -t ) 

21 p) " E y ^yaCC { t p )R hy (t m )R 2 , y (2t p ) [iZ) 

Here, the denominator can be regarded as an effective 
SLR function. This function takes into account that the 
contributions of the B and the M deuterons to the cos-cos 
correlation function are weighted with the initial corre- 
lations ciy C (t p ). Moreover, it includes that the relative 
contributions of the two species change due to diverse 
SSR, as is evident from the factors i?2,y(2i p ). 

In the Appendix, we discuss whether this approach en- 
ables correction for SLR. Our considerations show that, 
when two (or more) deuteron species contribute to the 
cos-cos correlation functions, elimination of relaxation ef- 
fects is successful only when three conditions are met: (i) 
Both species exhibit the same dynamical behavior, i.e., 
(f) B c (t m ; t p ) = (f>^f (t m ; t p ), which is fulfilled to a good ap- 
proximation in our case, see Sec. HVBl (ii) Appropriate 
evolution times are used so that z B c = z^f ~0. (iii) The 
data are corrected utilizing the effective SLR function de- 
fined in Eq. 112fl . The experimental determination of this 
function is described in the Appendix. In the following, 
we will mostly focus on evolution times, for which cor- 
rection for SLR is possible, i.e., (t m ; t p )~ F!^ (t m ;t p ). 
In these cases, the corrected data will be referred to as 
F§ c ', while we will use F 2 * to denote situations, when 
elimination of SLR effects is not successful. 

C. 7 Li NMR stimulated-echo spectroscopy 

Since the methodologies of 2 H and 7 Li NMR SE spec- 
troscopies are very similar, we restrict ourselves to point 
out some differences. In 7 Li NMR, the use of appro- 
priate three-pulse sequences allows one to measure F% , 
whereas F§ G is not accessible for 1 — 3/ 2. 56 Moreover, 
the shape of the quadrupolar coupling tensor is well de- 
fined in 2 H NMR, while it is subject to a broad distribu- 
tion in 7 Li NMR, leading to a damping of the oscillatory 
behavior of the initial and final correlations, in particu- 
lar, z ss = for sufficiently long evolution times. Further- 
more, in our case, 7 Li SLR is slow so that correction of 
the 7 Li NMR correlation functions for relaxation effects 
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is not necessary. Finally, we again emphasize that both 
ion and polymer dynamics can contribute to the decays 
of 7 Li NMR correlation functions and, hence, their inter- 
pretation is not as straightforward as for the 2 H NMR 
analogs. 



III. EXPERIMENT 



In addition to neat PPO, we study the polymer elec- 
trolytes 30:1 PPO-LiC10 4 , 15:1 PPO-LiC10 4 and 6:1 
PPO-LiC104. In the following, we will refer to these sam- 
ples as PPO, 30:1 PPO, 15:1 PPO and 6:1 PPO. In all 
samples, we use a deuterated polymer [CD 2 -CD(CD 3 )- 
0] n with a molecular weight M w = 6510g/mol. Strictly 
speaking, the polymer is poly (propylene glycol), having 
terminal OH groups. The sample preparation was de- 
scribed in our previous work4^ In DSC experiments^ 5 - 
we found T g = 203 ±2K for both PPO and 30:1 PPO 
and T g = 272 ±2K for 6:1 PPO. For the intermediate 
15:1 composition, we observed a broad glass transition 
step ranging from about 205 K to 245 K. In the litera- 
ture, two glass transition temperatures T 3i i «208K and 
T 9! 2«240K were reported for 16:1 PPO-LiC10 4 ^ 

The 2 H NMR experiments were performed on Bruker 
DSX 400 and DSX 500 spectrometers working at Lar- 
mor frequencies ujq/2tt of 61.4 MHz and 76.8 MHz, re- 
spectively. Two Bruker probes were used to apply the 
radio frequency pulses, resulting in 90° pulse lengths of 
2.2 /is and 3.5 /is. Comparing results for different setups, 
we determined that the results depend neither on the 
magnetic field strength nor on the probe, see Sec. IIV Al 
The 7 Li NMR measurements were carried out on the 
DSX 500 spectrometer working at a Larmor frequency 
of 194.4 MHz. In these experiments, the 90° pulse length 
amounted to 2.5 /is. Further experimental details can be 
found in our previous work4£ 

Appropriate three-pulse sequences P\ - t p - P 2 - t m - 
Pz - t are sufficient to measure sin-sin and cos-cos cor- 
relation functions for sufficiently large evolution times. 
However, for short evolution times, the stimulated echo, 
which occurs at t = t p , vanishes in the dead time of the 
receiver, following the excitation pulsesii 5 - This problem 
can be overcome, when a fourth pulse is added to re- 
focus the stimulated echo outside the dead time, lead- 



ing to four-pulse sequences Pi - t p - P2 - t m - P3 - A 
- We use three-pulse and four-pulse sequences to 

measure correlation functions for large (t p = 60 fis) and 
small (t p = 2 (is) evolution times, respectively. When we 
study the dependence of the cos-cos correlation functions 
on the evolution time, we apply a four-pulse sequence for 
all values of t p . In all experiments, we use suitable phase 
cycles to eliminate unwanted signal contributions. 57 
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FIG. 2: Correlation functions F% c (t m \t v = 60 fis) for neat 
PPO at various temperatures. The dashed lines are fits to 
Eq. with /3 = 0.33 and C = 0.01. 
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FIG. 3: Correlation functions F^ (t m \t v = 60 (is) for 6:1 
PPO-LiC104 at various temperatures. The dashed lines are 
fits to Eq. JTSJ with ,3 = 0.24 and C = 0.01. 



IV. RESULTS 

A. Temperature dependent 2 H NMR correlation 
functions 

In Figs. and we show temperature dependent cor- 
relation functions F§ (t m ; t p = 60 /is) for PPO and 6:1 
PPO. Due to the finite pulse lengths this value of t p cor- 
responds to an effective evolution time of about 62 /is, 
see Sec. 



IIV Bl so that the residual correlation for both 
deuteron species approximately vanishes and, hence, cor- 
rection for SLR is possible. Furthermore, we will show in 
Sec. IIV Bl that, for this large evolution time, F 2 (t m ]t p ) 
approximates the jump correlation function P 2 J (t m ), see 
Eq. ©. In Figs.|2and|3 it is evident that both materials 
exhibit nonexponential correlation functions. Moreover, 
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FIG. 4: Scaled cos-cos correlation functions for PPO, 30:1 
PPO-LiC10 4 , 15:1 PPO-LiC10 4 and 6:1 PPO-LiC10 4 (t p = 
60 jits). The temperatures are indicated. They are chosen so as 
to compare correlation functions with similar time constants 
r = 3.8— 5.8 ms. For the sake of comparison, the data was 
scaled according to Sf c (t m ) = F2 C (t m )/F2 C (t ) where t = 
5.6 ms w t . The dashed lines are interpolations with a function 
proportional to Eq. 1131 . The used stretching parameters are 
/3 = 0.33 for PPO, /3 = 0.3 for 30:1 PPO-LiC10 4 , = 0.2 for 
15:1 PPO-LiC10 4 and /3 = 0.24 for 6:1 PPO-LiC10 4 . 



for the 6:1 polymer electrolyte, the loss of correlation oc- 
curs at much higher temperatures in the time window of 
the experiment, indicating that the presence of salt slows 
down the polymer segmental motion. 

To study how the presence of salt affects the nonexpo- 
nentiality of the polymer dynamics, we compare correla- 
tion functions F^ ' (t m ;t p — 60 fis) for all studied mate- 
rials in Fig. 0] The respective temperatures were chosen 
so as to consider correlation functions with comparable 
time constants. For the sake of comparison, we actually 
show scaled correlation functions S% (tm)- We see that 
the nonexponentiality of the polymer dynamics is not 
a linear function of the salt content, but it is most pro- 
nounced for 15:1 PPO. Assuming that the rate of the seg- 
mental motion reflects the local salt concentration, these 
findings suggest that structural heterogeneities relating 
to the spatial salt distribution are most pronounced for 
the intermediate 15:1 composition, while the structure is 
more homogeneous for higher and lower salt contents. 

For a quantitative analysis of the correlation functions 
F2 (tm> t p — 60 /is), we fit the data for all compositions 
and temperatures with a modified KWW function 



(1 - C) exp 



C 



(13) 



At high temperatures, when the decay to the final state 
value can be observed in our time window, we obtain 
small residual correlations C = 0.01 ±0.01, in harmony 
with the theoretical value for this evolution time, see 
Fig. n Therefore, we keep this value fixed in the fits for 
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FIG. 5: Mean correlation times for (a) PPO, (c) 30:1 PPO- 
LiC10 4 and (d) 6:1 PPO-LiC10 4 . The time constants (rc> 
and (tj) result from interpolations of F§ S (t m ; tp — 2 /is) and 
Fa (t m ;t p =60/is) with Eq. JT3J, respectively. While the 
time constants shown as open squares and crosses were ob- 
tained for u)o/(2n) = 61.4 MHz, that displayed as solid squares 
and stars were determined for luo/(2-k) = 76.8 MHz. For 30:1 
PPO-LiC10 4 , the solid triangle is a data point resulting from 
experiments using a four-pulse rather than a three-pulse se- 
quence. In panel (b), we compare our results for PPO with 
correlation times r a from dielectric spectroscopy (I: Ref. IHsL 
II: Ref. The correlation times are shown on a reduced 

Tg/T scale to remove effects due to different glass transition 
temperatures, T g = 203 K and T g = 198 Kp^S The latter dif- 
ference reflects the fact that polymers with a smaller molecu- 
lar weight M m = 4000 g/ mol were investigated in the previous 
approaches. 25 ' 59 



the lower temperatures. For neither of the materials, we 
find evidence for a systematic temperature dependence of 
the stretching parameter (3. Specifically, we observe val- 
ues /3 = 0.33±0.02 for PPO, /3 = 0.30±0.02 for 30:1 PPO 
and /3 = 0.24±0.02 for 6:1 PPO. For the 15:1 composi- 
tion, determination of the stretching parameter is some- 
what ambiguous due to the extreme nonexponentiality. 
However, for all studied temperatures, a reasonable in- 
terpolation of the data can be obtained with /3 = 0.2, as 
will be discussed in more detail in Sec. IIV CI Hence, the 
values of the stretching parameters confirm that the non- 
exponentiality of the segmental motion depends on the 
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salt concentration. 

The mean correlation times can be obtained accord- 
ing to (r) = (t//3) T(l//3) from the fit parameters, where 
F(x) is the r function. In the following, we denote the 
mean time constants resulting for t p = 60 (is as mean jump 
correlation times (tj), since correlation functions for this 
relatively large evolution time approximate F^itm), see 
Sec. IIVBI The mean jump correlation times for PPO, 
30:1 PPO and 6:1 PPO are shown in Fig.0 We see that, 
for T g < T< 1.1 T g , the temperature dependence of (tj) is 
well described by an Arrhenius law. While we find similar 
activation energies i? Q = 2.7±0.2eV and £' a = 2.4±0.2eV 
for PPO and 30:1 PPO, respectively, there is a higher 
value E a = 4.1 ±0.2 eV for 6:1 PPO. The higher activa- 
tion energy for the latter composition is consistent with 
the slowing down of the segmental motion for high salt 
concentrations. In addition, we see in Fig.[5]that the cor- 
relation times are independent of the experimental setup, 
as one would expect. 

To obtain an estimate of the "NMR glass transition 
temperature" , we extrapolate the Arrhenius temperature 
dependence of the mean jump correlation time and de- 
termine the temperatures at which (r,/) = 100s. We find 
temperatures 205 K for PPO (T g = 203K), 208 K for 30:1 
PPO (T g = 203 K) and 283 K for 6:1 PPO (T 3 = 272K), 
which are in reasonable agreement with the calorimetri- 
cally determined values of T g . Thus, the slowing down of 
the segmental motion upon addition of salt reflects the 
shift of the glass transition temperature. 

Next, we focus on F 2 ss (i m ; t p = 2 /us) « F 2 (t m ). Since 
correction for SLR is not possible for sin-sin correlation 
functions, see Sec. Ill Bl this analysis is restricted to rel- 
atively high temperatures at which segmental motion is 
faster than SLR. In Fig. [5] we include the mean correla- 
tion times (tc), as obtained from fits ofFi ls (t m ;t p = 2iis) 
to Eq. I|13|) . We see that (tc) and (tj) exhibit a com- 
parable temperature dependence. For all compositions, 
(tc) is about one order of magnitude longer than (tj). 
Thus, a similar ratio (tj) / (tq) ~ 0.1 implies that the 
presence of salt does not affect the jump angles of the 
segmental motion. Using this ratio in Eq. J7|), we obtain 
a mean jump angle 7 = 15°. Hence, the segmental motion 
involves rotational jumps about small angles both in the 
absence and in the presence of salt. In Sec. UVBl the de- 
pendence of the mean correlation time on the evolution 
time, (r(tp)), will be studied in more detail. 

In dielectric spectroscopy, it was found that a Vogel- 
Fulcher behavior describes the temperature dependence 
of the a process of PPO in the presence and in the ab- 
sence of saltiS^SiS Thus, on first glance, it is surpris- 
ing that the correlation times follow an Arrhenius law in 
our case. For PPO, we show (tj(T)) and (t c (T)) to- 
gether with dielectric data r Q (T)22* 59 in Fig.^b). The 
time constants are plotted on a reduced T g /T scale to re- 
move effects resulting from the somewhat different glass 
transition temperatures of the studied polymers. Consid- 
ering the error resulting from this scaling, we conclude 
that our results are in reasonable agreement with the 
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FIG. 6: Correlation functions F$ (t m ;t p = 2 /is) for 30:1 
PPO-LiC104 at T = 226 K. The measurements were per- 
formed for complete recovery (£r = Is Ti,s) and partial 
recovery (tn = 3ms <Ti,s) of the magnetization after satu- 
ration. 



literature data. In particular, all the curves exhibit, if 
at all, a very small curvature in the temperature range 
0.9 < Tg/T < 1, rationalizing that an Arrhenius law en- 
ables a good description of the temperature dependence 
in this regime. At higher temperatures, the dielectric 
data show a pronounced curvature in the Arrhenius rep- 
resentation, as described by the Vogel-Fulcher behavior. 
In our case, the observation of unphysical prc-exponential 
factors To = 10~ 56 — 10~ 71 s implies a breakdown of the 
Arrhenius law at higher temperatures. Consistently, a 
previous SE study on the a process of a polymer melt 
reported prefactors tq « 10~ 77 s from Arrhenius inter- 
polations of temperature dependent correlation times>i& 
Hence, SE spectroscopy yields apparent activation en- 
ergies E a , characterizing the strong temperature depen- 
dence of the segmental motion close to T g , while devi- 
ations from an Arrhenius law become evident when fol- 
lowing T a (T) to higher temperatures. 



B. 



2 H NMR correlation functions for various 
relaxation delays and evolution times 



Next, we study whether different structural units of the 
PPO molecule show distinguishable reorientation behav- 
ior. Exploiting that T 1)M -CTi )B for PPO and 30:1 PPO 
at TK,T g r^ we record correlation functions for different 
relaxation delays tjt between the saturation of the mag- 
netization and the start of the SE pulse sequence. Then, 
the contribution of the M deuterons can be singled out in 
partially relaxed experiments, i.e., for tjjwTi^ <£.T± t B, 
while, of course, both deuteron species will contribute to 
the correlation functions if we wait until the recovery of 
the magnetization is complete. 

In Fig. El we display completely relaxed and partially 
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FIG. 7: Residual correlations from interpolations of cos-cos 
correlation functions for PPO at T = 222 K with Eq. JT3J. 
The top and bottom panel show data obtained for partial 
and complete recovery of the magnetization, respectively. The 
lines are theoretical values zf c /a,y C for 5m = 27t x39kHz (top 
panel) and <$g = 27r x 119 kHz (bottom panel) 



relaxed correlation functions for 30:1 PPO at T=226K. 
Evidently, s (i m ; t p — 2 fjs) is independent of the relax- 
ation delay tn, indicating that the "relevant" rotational 
jumps of the B and the M deuterons occur on essentially 
the same time scale. Here, the term "relevant" is used 
to denote the situation that the threefold jumps of the 
methyl group are not probed in the present approach. 
Performing an analogous analysis for neat PPO, we ar- 
rive at the same conclusion. Partially and completely 
relaxed correlation functions for larger values of t p will 
be compared below. 

In previous work on molecular and polymeric glass 
formers, it was demonstrated that more detailed in- 
sights into the reorientational mechanisms are possi- 
ble when F^ c (t m ' 1 t p ) is recorded for many evolution 
^i mes 1 5 4 i5 1 i6 J i7 1 i8 i p or example, information about the 
shape of the distribution of jump angles is available from 
the dependence of the mean correlation time on the evo- 
lution time, (r(t p )). Hence, it is desirable to apply 
this approach also to polymer electrolytes PPO-LiC104. 
However, in our case, the presence of two deuteron 
species hampers straightforward analysis of F^ c (t m ',t p ) 
for an arbitrary value of t p , see Sec. Ill Bl In the follow- 
ing, we test three approaches aiming to extract (r(t p )) in 
such situation. 

In approach (A), we measure partially relaxed corre- 
lation functions F§ (t m ',t p )- Then, the signal solely re- 
sults from the M deuterons and correction for SLR is 
possible for all values of t p when the data is divided by 
the relaxation function -Ri,Af(t m ), which can be obtained 
from SLR measurements. Thus, it is possible to deter- 
mine (r(t p )) from fits of the corrected data to Eq. (|13|1 . 
To demonstrate that this approach indeed allows us to 
single out the contribution of the M deuterons, we show 



experimental and theoretical values of the residual cor- 
relation F^ c (oo;t p ) as a function of t p in Fig. [7| The 
experimental values result from interpolations of the par- 
tially relaxed correlation functions F^ c (t m ;t p ) for PPO 
at T = 222 K. The theoretical values were obtained by 
calculating z%f (t p ) / a^f (t p ) for 5 M = 2n x 39kHzJ£ 
Obviously, the experimental and the theoretical curves 
F§ (oo;t p ) nicely agree, when we shift the measured 
data by 1.5 /is to longer evolution times. This shift re- 
flects the fact that, due to finite pulse lengths, the effec- 
tive evolution times are 1—2 fj,s longer than the set values, 
as was reported in previous work^i The good agreement 
of experimental and theoretical residual correlations in- 
dicates that approach (A) allows us to study the reorien- 
tational mechanisms of the methyl group axes. However, 
this approach is restricted to PPO and 30:1 PPO since 
T\^m and T^b differ less for higher salt concentrations, 45 
hampering a suppression of contributions from the B 
deuterons in partially relaxed experiments. 

In approaches (B) and (C), we focus on completely 
relaxed correlation functions. Then, we can include 6:1 
PPO, for which the salt should have the strongest ef- 
fect on the segmental motion. In both approaches, the 
measured decays for different evolution times are divided 
by the corresponding effective SLR functions, which are 
determined in concomitant SLR measurements, as is de- 
scribed in the Appendix. To emphasize that correction 
for SLR is not successful when Zg C (t p ), z^f (t p ) ^ 0, we 
denote the resulting data as F 2 *(t m ; t p ), see Eq. I|12f> . 

In approach (B), we extract (T(t p )) from interpolations 
of F2(t m ;t p ) with Eq. (|13|) . In doing so, we disregard 
that the residual correlation depends on the mixing time 
at t m rsTi m> see Appendix. In the limit t m —KX>, a resid- 
ual correlation F 2 *(oo; t p ) — Zg C (tp)/a^ c (tp) is expected 
based on Eq. 114|) . Figure[7]shows the residual correlation 
obtained from completely relaxed correlation functions 
F£(t m ;t p ) for PPO at T = 222K together with the the- 
oretical curve Zg C (t p ) / ag C (t p ) for S B = 2n x 119kHzi£i 
Shifting the measured points again by 1.5 /us, a good 
agreement of the experimental and the theoretical data 
is obtained, confirming our expectations for the residual 
correlation. Therefore, we use the theoretical values of 
F^oo; t p ) in a second round of fitting to reduce the num- 
ber of free parameters and the scattering of the data. 

In approach (C), we employ Eq. I|12l) to fit F£ (t m ;t p ). 
In doing so, we use the theoretical values of Zy (t p ) 
and a,y (t p ) together with the relaxation functions Ri, y 
and R,2.y from independent SLR measurements, see Ap- 
pendix. Considering our result that the relevant rota- 
tional jumps of the B and the M deuterons occur on 
very similar time scales, see Fig. |H| we utilize the same 
correlation function 4>(t m ;t p ) for both deuteron species. 
Assuming a KWW function, we have two free fit param- 
eters, r and j3. In Fig. [SI we show examples of Fjf (t m ; tp) 
for PPO at T = 222 K. Evidently, approach (C) enables 
a good description of the data for all values of t p . 

In Fig. El we present mean correlation times {r(t p )) 
resulting from the described approaches for PPO, 30:1 
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FIG. 9: Scaled mean correlation times (T(f p ))/(rc) for PPO 
at T = 222 K and 6:1 PPO-LiC10 4 at T = 297K and T = 301K. 
Results from approach (B) and approach (C) are included. 
For PPO, we compare data from two independent measure- 
ments (C1,C2). In addition, we show the scaled mean cor- 
relation times from partially relaxed experiments for PPO at 
T = 222K and for 30:1 PPO-LiC10 4 at T = 226K (t fl = 3ms), 
which were obtained using approach (A). See text for details. 



PPO and 6:1 PPO. Strictly speaking, the scaled data 
f(t p ) = (r(f p )) / (tc) are shown for the sake of compari- 
son. Inspecting the results of the partially relaxed experi- 
ments, it is evident that the rotational jumps in PPO and 
30:1 PPO are characterized by very similar curves T(t p ). 
Hence, the threefold methyl group axes show essentially 
the same reorientational mechanism in both materials. 
In particular, for both compositions, the dependence on 
the evolution time is weak for t p > 60 /is, indicating that 
the jump correlation function F% (t m ) is approximately 
measured. Thus, r{t p > 60 /is) provides an estimate for 
(t.j)/(tc)- Using a typical value t = 0.15, Eq. (JJJ yields 
a mean jump angle 7 = 18° for the rotational jumps of 
the methyl-group axes. 

Next, we focus on the outcome of the completely re- 
laxed experiments in Fig. |5J Although approaches (B) 
and (C) give somewhat different results, the overall shape 
of the respective curves f{t p ) is still comparable. In par- 
ticular, the findings for sufficiently large evolution times 
are in good agreement due to small values of the residual 
correlations Zy(t p ), see Fig. ^ We conclude that semi- 
quantitative information about the reorientational mech- 
anism is available in our case of two deuteron species. 
While the data allows us to compare the results for dif- 
ferent salt concentrations and to extract mean jump an- 
gles, we refrain from determination of the distribution of 
jump angles from further analysis of r(i p ), since the error 
is too big to resolve details of the shape of this distribu- 
tion. We estimate the data to be correct within a factor 
of two. The rough agreement of the curves r{t p ) for PPO 
and 6:1 PPO nevertheless indicates that the presence of 
salt, if at all, has a weak effect on the reorientational 
mechanism of the polymer segments, in harmony with 
the results in Fig. [S] In Fig. [5| we further see that r(t p ) 
weakly depends on t p for i p >40/is, justifying our argu- 
ment F§ c (t m ; t p = 60 /is) w F/ (t m ) . For both materials, 
we find r{t p > 40 /is) =0.1 — 0.2, corresponding to mean 
jump angles 7 = 15 — 21°. 

In Fig. 1101 we replot the results from the partially and 
the completely relaxed experiments for PPO at T — 222 K 
on a reduced St p scale. In this way, we can alleviate ef- 
fects due to the difference of the anisotropy parameters 
Smj^Sb- The good agreement of the completely and the 
partially relaxed data indicates that the relevant rota- 
tional jumps of the B and the M deuterons are similar 
not only with respect to the time scale, see Fig. but 
also with respect to the reorientational mechanism. Some 
deviations are expected since M and B deuterons con- 
tribute to the completely relaxed data and, hence, com- 
plete elimination of effects due to the different anisotropy 
parameters is not possible. 



C. 2 H NMR and 7 Li NMR correlation functions 
for 15:1 PPO-LiC10 4 

In Fig. we compare 2 H and 7 Li NMR correlation 
functions for 15:1 PPO. It is evident that the correlation 
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FIG. 11: Correlation functions for 15:1 PPO-LiC10 4 at the in- 
dicated temperatures: (a) S (tm', t p = 60 fts) from 7 Li NMR 
and (b) F^ c (t m ; t p = 60/j,s) from 2 H NMR. The dashed lines 
are fits with a KWW function, see Eq. ltT3l . where C = 0.00 
and C = 0.01 for 7 Li and 2 H NMR, respectively. The dotted 
lines are interpolations of the 7 Li NMR data for T — 253 K 
and T = 263 K with a two-step decay, (c) Mean correlation 
times (tj) determined based on the KWW fits of the 2 H and 
7 Li NMR results. We also include time constants obtained 
from measurements of Fjf s (t m ; i p = 60//s) in 2 H NMR. 



tions of about 15:1 PPO-LiC104) 25 i 3 ° one may speculate 
that a fit with two Kohlrausch functions gives a better de- 
scription of the 2 H NMR correlation functions. However, 
the limited time window does not allow us an unambigu- 
ous decision. In any case, the shape of the 2 H NMR 
correlation function shows that if fast and slow polymer 
dynamics exist the corresponding rate distributions will 
overlap in large parts so that there is no true bimodality. 
The 7 Li NMR correlation functions F^ s (t m ;t p = 60 /is) 
for T < 250 K are well described by a KWW function 
with /3«0.24. At higher temperatures, two-step decays 
are evident, suggesting the presence of two distinguish- 
able lithium species. However, since the interpretation of 
7 Li NMR data for polymer electrolytes has not yet been 
established, this tentative speculation deserves a careful 
check in future work, see Sec. Ill Al 

The mean time constants resulting from the KWW fits 
of the 2 H and 7 Li NMR correlation functions are com- 
piled in Fig. Hir e). At variance with the results for the 
other compositions, the mean jump correlation times (tj) 
extracted from the 2 H NMR data for 15:1 PPO do not 
follow an Arrhcnius law. For T>245K, the temperature 
dependence is described by _E a = 2.4eV, which is consis- 
tent with the apparent activation energies for PPO and 
30:1 PPO, see Fig.[SJ However, the variation of (tj) with 
temperature is weaker at lower temperatures. There are 
two possible explanations for the latter behavior. On the 
one hand, it is possible that a description using a single 
KWW function is not appropriate due to an existence 
of fast and slow polymer dynamics. On the other hand, 
it was demonstrated that the decay times of correlation 
functions do no longer reflect the shift of a distribution of 
correlation times when this distribution is much broader 
than the experimental time window as one may ex- 
pect for 15:1 PP04^ Rather, due to interference with 
the experimental time window, the decay times feign a 
weaker temperature dependence^ as is observed in our 
case. The mean correlation times resulting from the 7 Li 
NMR experiments are comparable to the corresponding 
2 H NMR data. This finding suggests a coupling of ion 
and polymer dynamics, as was inferred from visual in- 
spection of the correlation functions. However, we again 
emphasize that, in studies on the ion dynamics in poly- 
mer electrolytes, the interpretation of 7 Li NMR experi- 
ments is not yet clear. 



SUMMARY AND CONCLUSION 



functions for both nuclei decay on similar time scales, 
suggesting a coupling of ion and polymer dynamics. 
In 2 H NMR, reasonable interpolations of F^ c (t m ;t p — 
60 /js) are obtained for all studied temperatures when 
we use a KWW function with j3 = 0.2. However, closer 
inspection of the data reveals some deviations from a 
KWW behavior. Considering that a coexistence of fast 
and slow polymer dynamics was reported for composi- 



We demonstrated that 2 H NMR SE spectroscopy pro- 
vides straightforward access to both the time scale and 
the geometry of slow polymer segmental motion in poly- 
mer electrolytes. Here, we exploited these capabilities for 
an investigation of PPO, 30:1 PPO-LiC10 4 , 15:1 PPO- 
LiC104 and 6:1 PPO-LiC104. Comparison of the results 
for the different salt concentrations gives clear evidence 
that addition of salt leads to a strong slowing down of the 
polymer dynamics, while the geometry of the segmental 
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motion is hardly affected. 

The dependence of the 2 H NMR correlation functions 
on the evolution time yields quantitative information 
about the geometry of the segmental reorientationi^S 
This possibility results from the fact that the evolution 
time has the meaning of a geometrical filter. The jump 
angles of rotational jumps can be determined from the 
variation of the mean correlation time with the evolution 
time, (rftp)). 15 i 16 i 17 i 18 For PPO and 6:1 PPO, we found 
comparable curves (r(i p )), indicating that the polymer 
segments exhibit similar jump angle distributions in the 
presence and in the absence of salt. Furthermore, we ex- 
ploited that the mean jump angle can be determined from 
the ratio of the jump correlation time (tj) and the over- 
all correlation time (tc), which can be extracted from 
the 2 H NMR correlation functions in the limits t p — > oo 
and t p -> 0, respectively. For PPO, 30:1 PPO and 6:1 
PPO, a ratio of about 0.1 corresponds to a mean jump 
angle of ca. 15°. In the case of PPO, this value is con- 
sistent with results for various neat polymersj 3 *2iiLi2iAiii 8 . 
suggesting that small angle rotational jumps are typical 
of polymer segmental motion close to T g . For the poly- 
mer electrolytes, the agreement of the mean jump angles 
confirms that the reorientational mechanism is essentially 
independent of the salt concentration. 

In general, trans-gauche conformational transitions are 
attributed to the a process, implying the existence of 
large angle jumps. Thus, it may be surprising that the 
present and previous NMR studies on polymer melts and 
polymer electrolytes find small jump angles of 10 — 40°. 
To rationalize the latter behavior, the cooperativity of 
the a process was taken into account in the literature^ 
Specifically, it was argued that concerted conformational 
transitions of neighboring polymer segments lead to jump 
angles that are not well-defined and smaller than the an- 
gles typical of conformational transitions. 

Comparing the correlation times for the different mix- 
tures PPO-LiC104, we quantified the slowing down of 
the segmental motion as a consequence of the presence of 
salt. For PPO, the mean jump correlation time amounts 
to (tj) =1s at T= 212 K, whereas, for the highest salt 
concentration, 6:1 PPO, it has this value at T = 292K. 
Considering the strong temperature dependence of poly- 
mer dynamics close to T g , the difference AT = 80 K shows 
that addition of salt results in an enormous slowing down 
of the polymer segmental motion. This finding is consis- 
tent with both pronounced effects on the glass transition 
temperature and results from previous work. 32,58 

In the studied temperature ranges of T g < T < 1.1 T g , 
the temperature dependence of the mean jump correla- 
tion time (tj) follows an Arrhenius law. While we ob- 
served apparent activation energies E a ps 2.5 eV for the 
neat polymer and the moderate salt concentrations, a 
higher value E a w4.1 eV characterizes the segmental mo- 
tion in 6:1 PPO. Likewise, the mean correlation time 
(tc), characterizing the rotational correlation function 
of the second Legendre polynomial, shows a substantially 
stronger temperature dependence for 6:1 PPO than for 



PPO. Hence, the slowing down of the polymer dynam- 
ics for high salt concentrations is also reflected in higher 
apparent activation energies close to T g . For all compo- 
sitions, the pre-exponcntial factors of the Arrhenius laws 
have unphysical values, tq = 10~ 56 — 10~ 71 s, indicating 
that there is no simple thermally activated process over 
a single barrier. Consistently, a Vogel-Fulcher behavior 
can be used to describe the temperature dependence of 
the a process in a broader temperature ran g e ^i2S4SiS 
The 2 H NMR correlation functions for mixtures PPO- 
LiC104 decay in a strongly nonexponcntial manner, 
where the stretching depends on the salt concentration. 
For the correlation function F^ (t m ; t p = 60 /is), we ob- 
tained stretching parameters (3 « 0.33 for PPO, j3 ~ 0.3 
for 30:1 PPO, (3^0.2 for 15:1 PPO and /?«0.24 for 6:1 
PPO, which show no systematic variation with temper- 
ature in the range T g < T < 1.1 T g . Thus, the nonexpo- 
nentiality of the segmental motion is most prominent for 
intermediate salt concentrations of about 15:1, while it 
is reduced for both low and high salt content. Likewise, 
previous work reported that polymer dynamics is more 
nonexponcntial in the presence than in the absence of 
salti 32 i 61 We conclude that the additional structural het- 
erogeneities imposed by the addition of salt have a strong 
effect on the structural relaxation of the polymer melt. 

Comparison of 2 H NMR correlation functions for par- 
tial and complete recovery of the longitudinal magnetiza- 
tion provided insights into the reorientation of different 
structural units. With the exception of methyl group 
rotation, we found no evidence for diverse dynamical be- 
haviors of different structural units. Strictly speaking, 
the C-D bonds of the ethylene groups and the threefold 
methyl group axes show comparable rotational jumps 
with respect to both time scale and jump angles. 

All these results of the present 2 H NMR SE study are 
consistent with findings of our previous 2 H NMR line- 
shape analysis for polymer electrolytes PPO-LiC104. 45 
Specifically, the 2 H NMR spectra gave evidence that the 
presence of salt slows down the segmental motion. In 
addition, the 2 H NMR spectra indicated that, in large 
parts, the pronounced nonexponentiality of the polymer 
dynamics results from the existence of a broad distri- 
bution of correlation times G(logr). Then, the stretch- 
ing parameters (3 are a measure of the decadic width 
A of this distribution. In detail, the width is given by 
A ps 1.14(0.937/3 + 0.06)^ Based on this expression, we 
obtain widths of A = 3.3 and A = 5.4 orders of magni- 
tude from the stretching parameters for PPO and 15:1 
PPO, respectively. The latter result is in nice agreement 
with findings from our 2 H NMR line-shape analysis for 
the 15:1 composition. 45 There, the spectra were described 
using a logarithmic Gaussian distribution G(logr) with 
a full width at half maximum of 5-6 orders of magni- 
tude. We conclude that both 2 H NMR techniques yield 
a coherent picture of the segmental motion in polymer 
electrolytes. However, they highlight different aspects of 
the dynamics and, hence, a combined approach is partic- 
ularly useful. While 2 H NMR SE spectroscopy is a pow- 
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erful tool to quantify the time scale and the geometry 
of the segmental motion, 2 H NMR spectra allow one to 
demonstrate the existence of dynamical heterogeneities 
provided the distribution G(logr) is sufficiently broad. 

Having established that presence of salt slows down 
the segmental motion, one expects that the shape of the 
distribution of correlation times governing the polymer 
dynamics reflects features of the salt distribution in poly- 
mer electrolytes PPO-LiC104. We observed that the dis- 
tribution G(logr) is broadest for the 15:1 composition, 
implying a strong variation in the local salt concentra- 
tion for this intermediate composition, whereas the spa- 
tial distribution of the salt is more homogeneous for both 
high and low salt content. In the literature, it was pro- 
posed that there is liquid-liquid phase separation into 
salt-rich and salt-depleted regions for intermediate salt 
concentrationsi2Si2§422iaa Provided the rate of the segmen- 
tal motion is a measure of the local salt concentration, 
the existence of well-defined microphases should lead to 
a bimodal distribution of correlation times for the poly- 
mer dynamics and, consequently, to a two-step decay of 
the 2 H NMR correlation functions. Although there are 
some indications for the presence of a two-step decay, a 
KWW function still provides a satisfactory interpolation 
of the experimental data for 15:1 PPO-LiC104. Hence, 
our results are at variance with a bimodal, discontinuous 
distribution of correlation times and, thus, of local salt 
concentrations unless the contributions from salt-rich and 
salt-depleted regions overlap in large parts, while they are 
consistent with the existence of a continuous G(log r). In 
particular, the 2 H NMR spectra and correlation functions 
for 15:1 PPO are in harmony with the picture that short 
ranged fluctuations of the local salt concentration lead 
to a large structural diversity and, thus, to a broad con- 
tinuous rate distribution for the segmental motion. For 
example, our results are consistent with a reverse Monte- 
Carlo model of 16:1 PPO-LiC104, which features salt-rich 
and salt-depleted regions on a short length scale of about 
10 AjjS corresponding to about two monomeric units. 

In addition, we performed 7 Li NMR SE experiments to 
investigate the lithium ionic motion in 15:1 PPO-LiC104. 
For polymer electrolytes, a strong coupling of ion and 
polymer dynamics was proposed in the literatures^ 
Consistently, we found nonexponential 7 Li NMR correla- 
tion functions, which agree roughly with the correspond- 
ing 2 H NMR data. Hence, our results are not at variance 
with the common picture of charge transport in polymer 
electrolytes. However, we emphasize that the present 
findings cannot confirm a coupling of ion and polymer 
dynamics beyond doubt, since the interpretation of 7 Li 
NMR studies on lithium ionic motion in polymer elec- 
trolytes needs to be further clarified. Specifically, it is 
not clear to which extent the time dependence of the 7 Li 
NMR resonance frequency reflects the dynamics of the 
ion and the rearrangement of neighboring polymer seg- 
ments, respectively. 

From the viewpoint of NMR methodology, we inves- 
tigated whether 2 H NMR SE spectroscopy provides well 



defined information when there are contributions from 
two or more deuteron species, exhibiting different relax- 
ation behaviors. We showed that problems can arise 
when at least one of the deuteron species exhibits fast 
SLR so that the decay of the 2 H NMR correlation func- 
tions due to relaxation interferes with that due to molec- 
ular dynamics. It was demonstrated that a successful 
correction for SLR requires that the different deuteron 
species show similar molecular dynamics. Moreover, it 
relies on the choice of appropriate evolution times t p . 
While the effects of SLR can be eliminated for evolution 
times t p , for which the residual correlations z(t p ) vanish 
for all deuteron species, correction for relaxation effects is 
not possible in the general case. Then, a more elaborate 
analysis can yield some insights, however the experimen- 
tal error increases. 



VI. APPENDIX 

Here, we discuss whether relaxation effects can be elim- 
inated when the experimental data are divided by an ef- 
fective SLR function, see Eq. l(T2")l, We assume that the 
same correlation function 4>(t m ;t p ) characterizes the rel- 
evant rotational jumps of the B and the M deuterons, as 
is suggested by experimental findings, see Fig. Then, 
short calculation shows that -F 2 *(i m ;i p ) = (j)(£ m \t p ) for 
Zfrf(t p ) = Zg C {t p ) = 0. Consequently, the cos-cos cor- 
relation functions provide straightforward access to the 
molecular reorientation when we use evolution times, for 
which the residual correlations of both deuteron species 
vanish. However, such correction for SLR is not possi- 
ble for an arbitrary value of t p . To illustrate the effects 
it suffices to consider the case i?2,s(2t p ) = i?2,Af (2£ p ). 
First, it is instructive to study the residual correlation of 
F2{t m ]t p ). For t m ^>r, we obtain from Eq. l(TT|l 

7?*(t \ — z M C fa)Pi(^™) + z B C (tp) / 1zL \ 

2( ™ ' P) " aCf(t p ) Pl (t m ) + a cc(t p ) (Mj 

where pi(t m ) =Ri,M(t m )/Ri,B(tm)- In our case Ti,m < 
Ti b, this function decays from p\{t m Ti m) = 1 to 
Pi(t m ^ Ti,A/) = and, thus, the residual correlation 
shows a crossover from [zm+zb]/[cim+<ib] to zb/cib at 
t m ~Ti Mi where we skipped the superscripts. When the 
decays due to molecular dynamics and SLR are well sep- 
arated, i.e., for t<cTijv/j this crossover is not relevant 
since the analysis of the correlation functions can be re- 
stricted to mixing times t m -^LTi.m- However, in the case 
T~2~i,Afi which is often found in the experiment due to 
broad distributions of correlation times, the crossover of 
the residual correlation at t m ~T\^m interferes with the 
decay due to molecular dynamics. Then, our approach 
does not enable an elimination of relaxation effects. 

To make the effects more clear, we show several calcu- 
lated correlation functions in Fig.^| For the calculation, 
we assume that the molecular rotational jumps are de- 
scribed by a KWW function 4>(t m ) =exp[— (i m /3ms) - 45 ]. 
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FIG. 12: Calculated cos-cos correlation functions for t p 
I U fis and t p — 65 /j,s. For the calculations, we utilized Eqs. @ 
and JTT), where <j> y (t m ; t p ) =exp[— (t m /3 ms) 0,45 ], R,2, v (2t p ) = 
1, Ri,M(t m ) = exp( — t m /10ms) and i?i,s(t m ) = exp(-i m /3 s). 
The values of the initial and the final correlations were com- 
puted for Sb = 27T x 120 kHz and <5a/ = 27r x 40 kHz accord- 
ing to Eqs. @ and (1101 . We show the correlation functions 
(t m ',t p ) and £§ (tm',tp), differing due to the action of 
SLR, together with the normalized counterparts F§ c \t m ;t p ) 
and J-§ C (tm\ t p ). Moreover, we include the correlation func- 
tions F2 (i m ;i P ) calculated according to Eq. 11211 . 

Further, Sb = 27rxl20 kHz and 5m = 27rx40 kHz are used to 
compute the initial and the residual correlations. Utiliz- 
ing these data, E^ c is obtained from Eq. JHJ. To mimic 
the damping due to SLR, we calculate according 
to Eq. 1|11|> . In doing so, we suppose exponential relax- 
ation functions Ri t y(t m ) characterized by T\^m — 10ms 
and Ti b = 3 s. Thus, the decays due to dynamics and 
relaxation interfere for the M deuterons. In Fig. ^1 the 
effect of SLR is evident from comparison of E^ and 
£2° ■ Normalizing these data, we obtain F^ and T^ C ' ■ 
Finally, Eq. (|12(l is used to compute F 2 *. If our approach 
enables elimination of SLR effects, F 2 *(t m ) and (t m ) 
will be identical. 

In Fig. 1121 we compare the various correlation func- 
tions for the evolution times t p — 10 fxs and t p = 65 fj,s. For 
the latter evolution time, the condition Zg = z M c « is 
met, see Fig. ^ so that elimination of relaxation effects 
should be possible. This expectation is confirmed by the 
nice agreement of F^tm) and F 2 cc (t m ). By contrast, for 
t p = 10 fis, wc find z^f > Zg C > 0. Then, F^tm) clearly 



differs from Ff c (i m ), indicating that correction for SLR 
is not successful. In particular, F^itm) exhibits a residual 
correlation Zg C /a§ c , which deviates substantially from 
the value resulting in the absence of relaxation effects, 
see Eq. (fTIj) . 

For an experimental determination of the effective SLR 
function, we measure the SLR using the saturation recov- 
ery technique, where the longitudinal magnetization is 
read out with a solid-echo so that SSR reduces the height 
of the solid-echo. For an evolution time t p in the SE 
experiment, we perform a concomitant SLR experiment 
choosing the solid-echo delay t e such that the total time 
intervals during which transversal magnetization exists 
are equal in both measurement. In detail, we use t e ~t p 
or t e =t p + A, depending on whether a three-pulse or a 
four-pulse sequence is applied in the SE experiment, see 
Sec. IIIII Then, in both the SLR and the SE experiment, 
the contribution of each deuteron species is proportional 
to the " effective" magnetization m y R2^ y (2t e ). In particu- 
lar, the recovery of the longitudinal magnetization is the 
sum of the contributions m y R2^ y (2t e )Ri, y (t m ). Multipli- 
cation of these contributions with the respective theoret- 
ical values of a y c (t p ) yields the effective SLR function. 

At the studied temperatures, the recovery of the longi- 
tudinal magnetization occurs in two steps, which can be 
attributed to the B and the M deuterons, respectively^ 
Therefore, it is straightforward to determine the relax- 
ation functions Ri jV (t m ) and the relative contributions 
of the two deuteron species. When we study the depen- 
dence of F§ (t m ]t p ) on the evolution time and, hence, 
when we use various solid-echo delays t e = t p + A in the 
concomitant SLR measurements, we find that Ri t y(t m ) 
does not exhibit a systematic variation with t e , as one 
would expect. Therefore, we use the average values for 
the calculation of the effective SLR functions. However, 
the relative amplitudes of the two steps depend on the 
echo delay, indicating that the B and the M deuterons 
show different SSR behavior. For PPO at T = 222K, the 
contribution of the M deuterons increases from 55% for 
t e = 25/is to 65% for t e = 85/is. Hence, the M deuterons 
exhibit slower SSR, consistent with a smaller anisotropy 
parameter 6. When extending the evolution time in SE 
experiments, this discrepancy in SSR means that the rel- 
ative contribution of the B deuterons decreases with re- 
spect to the theoretical value, a|/(a|+a|, f ). 
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